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Discussion of 
“EXAMPLES OF TIMBER STRUCTURE FAILURES” 


by Michael N. Salgo 
(Proc. Paper 544) 


MICHAEL N. SALGO,! M. ASCE.—Dr. Stern, in his discussion, inquires as 
to the type of split rings, flat face, bevelled on inner or both faces, for the 
structures reported on. All of the structures reported on had the standard 
split ring in use at that time namely, the split ring with a bevelled inside and 
flat outside face. In view of this, the desired information on field experiences 
with flat sided versus bevelled split ring connectors is not available from the 
writer. 

Mr. Gloss, in his discussion, sets forth some of his experiences and raises 
certain questions. Comments on these will be set forth in subsequent para- 
graphs. 

Mr. Gloss apparently is of the opinion that the paper and parts of the con- 
clusion may be construed as placing the blame for failures encountered on 
truss types or fastening means, rather than on the errors of design, fabrica- 
tion or maintenance. This was definitely not the intent of the paper. Quite the 
contrary, the writer took particular care, in discussing the failures, to point 
out the respective causes which were basically design and/or construction 
and/or maintenance. Certain types of framing did have higher incidence of 
failures than others and these were pointed out. 

As Mr. Gloss indicates, the writer is of the opinion that the over-all per- 
formance record of timber trusses during World War II was outstanding and, 
as stated in the paper, “the availability and use of long-span timber structures 
during World War II is considered to have been a major contribution to the 
war effort.” However, failures and difficulties did occur and from these cer- 
tain lessons can be learned. 

One of the major sources of difficulties with parallel chord trusses during 
World War II was the fact that actual deflections were considerably in excess 
of that which could have been predicted on the basis of available design data 
for deformation of split ring joints under load. Subsequent design data, as 
promulgated, increased the expected deflections for timber connector type 
trusses. However, by this time, many timber truss structures had been built. 
Results of this manifested itself in such items as lack of clearance for bridge 
cranes, binding of large hangar doors, etc. These structures were cambered 
on the basis of the then available design data. Where bowstring trusses were 
used, the effect of the deformation data was not as critical as in the parallel 
chord trusses. 

Mr. Gloss indicates particular interest in the discussion of secondary 
stresses. Although recognizing the presence of the secondary stresses, he 
questions the validity of most computations to determine their magnitude, 
since such computations involve many critical assumptions. It is the writer’s 
opinion that the assumptions are not necessarily as critical as Mr. Gloss sets 


1. Principal Asst. to Dist. Civ. Eng., Third Naval Dist. Headquarters, New 
York, N. Y. 
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forth, since the same assumptions that satisfactorily predict the total deflec- 
tions (including creep effects) can also be used to predict member distress. 
Figure 19 illustrated the results of such an analysis, which were subsequently 
used as the basis for necessary repairs and strengthening. However, these 
analyses are often quite difficult, time consuming, and beyond capabilities of 
ordinary design procedures; for this reason, inclusion of certain limiting fac- 
tors in building codes are considered warranted. 

There appears to be an honest difference of professional opinion concern- 
ing the proposed restrictions for trusses other than bowstring where unsea- 
soned lumber is to be used. To further belabor this point at this time does 
not appear warranted. However, at the time this paper was formally present- 
ed at the Annual New York Convention of the ASCE on October 21, 1953, the 
inclusion of the proposed restrictions was still pending. A decision was sub- 
sequently made not to have the restrictions in the AITC-Timber Construction 
Standards. The writer, however, still considers it to be good practice and so 
recommends. 

In closing, the writer wishes to restate the basic purpose of the paper 
which was to present certain actual observed experiences, particularly with 
reference to data on failures of timber structures observed from which cer- 
tain lessons can be learned; and the final conclusion that “TIMBER STRUC- 
TURES PROPERLY DESIGNED, CONSTRUCTED, AND MAINTAINED HAVE 
TAKEN THEIR PLACE AS MAJOR COMPONENTS IN OUR PRESENT DAY 
CONSTRUCTION PROGRAM.” 


at 
| 
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Discussion of 
“AN INFLUENCE LINE ANALYSIS FOR SUSPENSION BRIDGES” 


by David J. Peery 
(Proc. Paper 558) 


DAVID J. PEERY,! M. ASCE.—The discussions include significant supple- 
mentary theory, as well as substantiation of the writer’s basic influence line 
procedures. The contributions of the discussers are appreciated. Most of 
the discussions treat the simpler case of the bridge with hinged stiffening 
trusses. The obvious questions concern the accuracy and the relative sim- 
plicity of the influence line procedure, as compared with the customary de- 
flection theory. 

Mr. Sih questions whether the influence lines are applicable to unstiffened 
bridges, and feels that they give a “misleading concept of the function of the 
stiffening truss.” The influence lines shown in Fig. 16 include the complete 
range of values from cQ = 0 (elastic theory) to c2 =-o(unstiffened bridge). 
These influence lines are readily visualized as deflection curves, whereas the 
conventional long equations of exponential or hyperbolic terms do not present 
a simple physical picture. 

Mr. Sih also questions the assumption of superposition. In reality, this 
assumption is identical with that used in all deflection theory analysis, and 
influence line ordinates can be calculated to any desired degree of precision, 
by successive corrections of the cQ values. The basic differential equation of 
the deflection theory, 7"' - c2n = -M/EI is non-linear, since c is not constant. 
It is not feasible to integrate this non-linear equation, so c is always assumed 
constant during integration, thus linearizing the differential equation. Super- 
position procedures apply for solution of linear differential equations, but not 
for non-linear problems. 

Dr. Karol presents simple equations for the variation of H, M, and V with 
cQ. When these equations are used, it is necessary to calculate influence 
lines for only one value of cQ. Values can then be corrected to correspond to 
the true value of cQ for each particular loading. The computation work is ob- 
viously reduced from that required to calculate two sets of influence lines. 

Mr. Paulbaum extends the influence line procedure to the calculation of 
shears in the stiffening truss. This is an important contribution, since two 
segments of the span may be loaded and the influence lines provide the only 
convenient method for finding the correct loaded length. 

Professor Franciosi substantiates the author’s procedure of finding in- 
fluence lines as deflection curves, or by using the theorem of Betti. He proves 
the validity of this procedure for suspension bridges and arches, provided HR 
remains constant. He also proposes a method of successive approximation, 
which is applicable to arches or suspension bridges. The successive approxi- 
mations would probably converge more slowly for suspension bridges than for 
arches because of the higher values of c2. The equations for M are also 
more sensitive to changes in cQ than are the equations for H. 


1. Prof. of Aeronautical Eng., Univ. of Michigan, Ann Arbor, Mich. 
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Several of the discussers requested more design data and further exam- 
ples of the computational procedure. Charts and tables will therefore be pre- 
sented for bridges with simply supported trusses. Most bridges of this type 
are somewhat similar, with most of the flexibility contributed by the main 
span. The side spans are always considered as unloaded or as uniformly 
loaded for the entire span; hence, influence lines are not required for the side 
spans. It will be assumed that 80 percent of the denominator D of the H- 
equation (Eqs. 12 and 13) is contributed by the main span term, which will be 
designated as Dp, and the remaining 20 percent of D is contributed by the side 


L 
span terms 2D. and the cable stretch term ix This assumption is stated 
c 
in the following equations, which also define Dp and Dg. 
L 
s 
D= D, + 2D, + AE, (62) 


D 1.25 D, (63) 
et 
Dy =~ 92 J (64) 


| h (65) 


The influence lines will be calculated for bridges with these proportions. The 
final bending moments will then be corrected to the true ratios of D/D,. 

The value of Dp (or Dg) is obtained by substituting 7), from Eq. 11 into 
Eq. 64, and integrating. 


2 
p 1 2 | 66) 


For conversion to exponentials, tanh cQ/2 = 1)/(e° 1). Values of 
DoEI/{2Q or DgEl;/f;* 24 are listed in Table 6 for various values of cQ. 
These terms are also multiplied by c2Q2 in order to permit easy interpolation 
at the higher values of cQ. 

The H-influence line is now obtained by substituting D=1.25 Dp, and values 
of 1}, from Eq. 11 intoEq. 13. The H-influence lines for 9 = 10f are plotted 
in Fig. 16 for various values of cQ. The H-influence lines are seen to have 
little variation in shape, since the maximum ordinate varies from 1.5625 for 
cQ = 0 to 1.5000 for cQ =~o, or only 4 percent. This confirms Dr. Karol’s 
statement that the H-influence lines are affinely related. For other ratios of 
{/Q and D/D,, the H ordinates should be multiplied by 2/10f x 1.25 D,/D. Dr. 
Karol’s Eq. 37 depends on a typical variation of the ratio D,/D withcQ. The 
values of H decrease with higher values of cQ, since the ratio D)/D decreases. 
This effect is not shown in Fig. 16, which shows only the relative shapes of 
the H-influence lines. 

For a uniformly distributed live-load p extending a distance k' from the 
left end of the span, the values of H/pQ are obtained from the areas under the 
H-influence lines. These values are shown in Table 7. For live loads over 
the right side from k' toQ, the values in Table 7 may be subtracted from 
1.000, since H/pQ = 1.000 for a full live load. Any partial live loading may be 
considered by superposition. 
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Values of m)/f may be calculated from Eq. 10. For convenient interpola- 
tion, these values are multiplied by c2Q2 and listed in Table 8. The values of 
Mp obtained from Eqs. 8 and 9 depend only oncQ. The values of mp vary 
directly with the sag f, and the values of H vary inversely withf. The product 
Hm», is therefore independent of f, and varies only with cQ, when D = 1.25 Do: 
Consequently, the bending moment influence lines may be calculated from 

Eq. 14 as functions of cQ only, if D is assumed equal to 1.25 Dp. These bend- 
ing moment influence lines are plotted in Figs. 16 to 18. For normal ratios 

of a , the influence lines of Figs. 16 to 18 are sufficiently accurate to deter- 
mine the final load position for maximum positive or negative bending moment. 

The positive and negative areas under the influence lines of Figs. 16 to 18 
may be calculated by integration of Eqs. 8, 9, and 11. These values are plotted 
as cQx M/pQ? in Figs. 19 and 20. The factor c Q is included so that the 
curves may be plotted to a larger scale. 

The bending moments calculated for D/D, = 1.25 must now be corrected to 
the true ratio of D/Dg. If the value from Table 7 is designated as Hp, the in- 
crement §H produces a bending moment correction of 6H x mp. 
1.25 D, Q 

-1)— 
10f 


This value of 6H may be added to Ht to correct the positive main span bend- 
ing moments. 


(67) 


M = (Ht + 6 H)m), (68) 


The maximum negative bending moments occur when the side spans are fully 
loaded. A full live load in one side span produces a cable component which is 
designated by H,. The negative bending moments in the main span have the 

following correction: 


M = (2H, + H, + 6H)m, (69) 


The cable comp:-nent for full live load in one side span is obtained by integrat- 
ing Eq. 13 over one side span and substituting D, from Eq. 65. The value 
Hy = w1232/s is also substituted. 


P D, 
H, = Hy, (70) 


1 


The side span bending moments, Mg, are calculated directly from Eq. 10. 
For negative bending moments in one side span, the value of H is calculated 
for the other spans loaded at minimum temperature. 


M = Hm), (71) 


For positive bending moment in one side span, only that span is loaded at 
maximum temperature. The uniform load p produces a moment proportional 


2 
to the moment my by the ratio - PQ, or - Hy = The total positive bend- 
8f; 


ing moment is 


M = ( Hy w, * 4s H,) m, (72) 


The values of m» can be interpolated from Table 8. 
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Numerical Application 


The following structure will be analyzed. The dimensions are 2 = 3800 ft, 
Q, = 1800 ft, f = 350 ft, fy = 78.53 ft, L, = 9140 ft and Ly = 8482 ft. The flexi- 
bility properties for one side are EI = Ely = 1884 x 106 ft2 kips, and Ls/A Eg = 
0.000880 ft/kip. The load per cable is w = wy = 4.7 kip/ft and p = 1.0 kip/ft. 
The values t = +t 60°F, w = 0.0000065 and Hy = wQ2/ef = 24,240 kips are used. 

The main span bending moments will be interpolated between calculated 
values for cQ = 14 and cQ = 15,but only the calculations for c 2 = 14 are shown. 
For cQ = 14, cyQq = 14 x 1800/3800 = 6.63, and H, = c?EI = 25,570K, From 
Table 6, Do = 5.053 {*/H, Q = 0.006370, and Dg, = 4.314 x f12/H, 2, = 0.000578. 
From Eq. 62, D = 0.008406 ft/kip. The ratio D/D, is 1.320. For a tempera- 
ture change, Hy =  wtl,/D = * 393 kips. For full live load in the main span, 

D 
K s 
. For full live load in one side span, H, = Hy, os ie 


D 
3908 
354K. Maximum H = 3908 + 2 x 354 + 393 = 5009 kips. 

The bending moments are calculated in Table 9. The load positions are ob- 
tained from the influence lines. The main-span bending moments are first 
obtained from Figs. 19 and 20, which represent the positive and negative 
areas under the influence lines. These moments are designated as Mj in 
lines 5 and 11 of Table 9. The moments are then corrected for the D/Dp 
ratio, and for temperature and side span loads. The moments calculated in 
lines 7 and 13 are for c2 = 14. The moments shown in lines 8 and 14 are in- 
terpolated from similar calculations for cQ=15. The side span bending mo- 
ments are calculated for the correct values of ©42y in lines 16 and 18. 


SUMMARY 


The original paper presented methods for obtaining influence lines for 
suspension bridges with any number of spans, and with hinged or continuous 
stiffening trusses. The procedure permits a direct determination of the 
critical loading conditions, and permits consideration of any type of live load 
distribution, without the usual restriction of uniform loading. Influence lines 
are obtained by superimposing simple deflection curves, thus permitting a 
designer to detect gross errors. The customary long, abstract equations are 
avoided, Any desired degree of accuracy may be obtained by calculating in- 
fluence lines at small intervals of cQ, although engineering accuracy may al- 
ways be obtained with only two values of c 

For two-hinged suspension bridges, influence line procedures had previ- 
ously been used by Rode, 4 Selberg,° Asplund,® and Karol,? all of whom pre- 
pared tabular data. Similar charts and curves are readily prepared for con- 
tinuous bridges since the equations (18 and 19) are simple. However, tables 
covering the complete range of values appear too voluminous for this paper. 

The influence lines constructed for D = 1.25 Dp appear adequate for most 
two-hinged bridges. For bridges of unusual proportions, the influence lines 
may be easily corrected to give the exact load position. Similar influence 
lines may be constructed for shears. 


Acknowledgment: The data used in the numerical application correspond 
to preliminary dimensions of the Mackinac Straits Bridge, and were provided 
by Dr. D. B. Steinman, M. ASCE. 

4, 5, 6, 7. Footnotes in Proc. Sep. 558. 
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Table 8 
Values of -c22 2m)/f 


0.1.2 0.2 2 0.39 0.42 


1.842 3.126 3.969 4.445 
2.195 3.672 4.611 5.133 
2.519 4.150 5.156 5.701 
2.817 4.567 5.609 6.161 
3.092 4.932 5.987 6.529 
3.346 5.250 6.300 6.824 
3.584 5.530 6.560 7.059 
3.808 5.780 6.778 7.245 
4.014 6.001 6.961 7.393 
4.400 6.373 7.244 7.606 
4.746 6.672 7.447 7.745 
5.056 6.916 7.594 7.831 
5.337 7.112 7.701 7.891 
5.590 7.272 7.779 7.928 
5.820 7.405 7.837 7.952 
6.024 7.512 7.880 7.968 
6.213 ».602 7.911 7.979 
6.385 7.674 7.934 7.986 
6.678 7.782 7.964 7.994 
6.920 7.854 7.980 8.000 
7.854 7.998 8.000 8.000 
8.000 8.000 8.000 8.000 


Tabular values may also be used to calculate "= (y - m,)/Hp 


0.5.2 
3.0 4.599 
3.5 5.302 
4.0 5.873 
4.5 6.331 
5.0 6.695 
5.5 6.981 
6.0 7.205 
| 6.5 7.380 
7.0 7.517 1 
8.0 7.706 
9.0 7.824 
10.0 7.892 
11.0 7.934 
12.0 7.960 
13.0 7.976 
14,0 7.984 
15.0 7.991 
16.0 7.994 
18.0 8.000 
20.0 8.000 
40.0 8.000 
oo 8.000 
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Discussion of 
“ULTIMATE SLOPES AND DEFLECTIONS” 
A Brief Limit for Design 


by George C. Ernst 
(Proc. Paper 583) 


GEORGE C. ERNST,* M. ASCE.—Mr. Diwan’s discussion is based in its 
entirety upon the rejection of the completely idealized relationship between 

M and W for conditions prior to collapse and a questioning of the value of the 
“unit rotation diagram.” 

With regard to the M and w relationship, the true curve is a function of 
steel ratio, local bond development and failure, steel properties in tension, 
and concrete properties in tension and compression. Up to initial tensile 
cracking not exceeding a unit deformation in the concrete of about 0.0002, a 
reinforced concrete beam does not materially differ from a plain concrete 
beam insofar as slope and deflection relationships are concerned. Subsequent 
to tensile cracking, the tensile properties of the concrete have a gradually 
decreasing effect until local failure of bond occurs between cracks, thereby 
resulting in steel deformations conforming to Eg = 30,000,000 psi over sub- 
stantial portions of the beam. Crack formation during the early stages of 
loading does not produce accumulated rotations as calculated with an E, = 
30,000,000 psi as Mr. Diwan has assumed. From the development of yield in 
the steel at M, to the ultimate moment M,,, the shape of the curve is domi- 
nated by the compression properties of the concrete and the steel ratio! for 
any given grade of steel. In general, progressive local bond failure in the 
vicinity of tensile cracks begins at a unit deformation of about 0.0002 in the 
concrete and a steel stress at the crack of about 15,000 psi. Figure 13 com- 
pares the M -wW curves for three different assumptions as to the effective Eg, 
the intermediate one being most representative of the true variation. Figure 
14 compares the slopes and deflections for a simply supported beam with con- 
centrated load at midspan, for which the curvature from M, to M, in Mr. 
Diwan’s Type 2 curve has been used with each of the three assumptions for ef- 
fective Es. Obviously, the assumption of an infinite E, merely expresses the 
upper limit of stiffness for the tensile face, and the value of 30,000,000 psi for 
E, expresses the lower limit up to yield. 

Mr. Diwan’s belief that the curvature from M, to M, will provide a realis- 
tic estimate of the length of the plastic portion in the vicinity of a “peaked” 
My for purposes of moment redistribution in a continuous structure can be 
shown to be quite unwarranted. Unfortunately, Mr. Diwan has selected rela- 
tively favorable conditions for redistribution of moments in his examples. 
However, taking the values of Mp and M, from his table, the difference in 
moment is about 7% of My for a drop of 87% in unit rotation for AB, and 6.4% 
for a 93% drop for BC and CD. In the vicinity of a linearly peaked My, there- 
fore, 7% (for AB) or 6.4% (for BC or CD) of the distance to zero moment on 
either side of M, would be affected by the unit rotations between M, and My. 


* Prof. and Chairman, Dept. of Civ. Eng., Univ. of Nebraska, Lincoln, Nebr. 
10. “Inelastic Behavior of Reinforced Concrete Members subjected to Short- 
Time Static Loads,” L. H. N. Lee, ASCE Proc.-Separate No. 286, 1953. 
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For Mr. Diwan’s loading 2 on his single bay frame, the required length over 
which v/, tow, must be accumulated to produce ¢¢- is 12.0 inches, whereas 
the length actually affected from Mp to My is only 6.4% of 120 inches or 7.7 
inches. Apparently Mr. Diwan is depending upon some other unstated property 
to provide for the total required angle change. Figure 15 illustrates the effect 
of varying the steel ratio and concrete strength upon the minimum length over 
which plastic rotation must occur at the supports of a fixed end beam under 
uniform load in order to assure complete moment redistribution. It should be 
recalled that yw markedly decreases when the steel ratio is increased or the 
concrete strength reduced, !! thereby reducing the total angle change available 
for moment redistribution by virtue of the portion of the beam affected by M, 
to My. Contrarily, this in turn increases the required minimum length of 
plastic hinging as is shown in Figure 15. 

All this adds up to the fact that it should be reasonably obvious that moment 
redistribution cannot be significantly dependent upon the difference between 
Mr. Diwan’s Type 1 and Type 2 curves in his Figure 8, but rather on the dis- 
tribution of yy, as controlled by the distance throughout which steel yield can 
occur due to length of bearing, width of column, shape of knee, etc. The effect 
of these factors can be determined only by laboratory tests of full-sized mem- 
bers and structures. 

As to the value of the ‘unit rotation diagram,” it is undoubtedly realized by 
all readers that the M/EI diagram is a unit rotation diagram based upon a 
linear relationship between stress and strain. The M/EI concept for those 
portions of a structure below first yield was not presented because, for the 
purpose of the paper, it was felt best to avoid the obscuring of the purely geo- 
metrical relationship between deformations, slopes, and deflections. The 
author does not believe that the unit rotation diagram, either in its basic form 
or converted to M/EI form, will have any significant usage in ultimate load or 
limit design. It has a rather limited usefulness even for the so-called “elas- 
tic” method of analysis. The problem is essentially that of obtaining a funda- 
mental understanding of structural response throughout the loading history of 
a structure, and a very necessary part of the problem is the experimental 
determination of the extent of distribution of unit rotations beyond first yield. 
The paper was intended to draw attention to this feature as well as the mo- 
ment-slope-deflection relationship, and the author is indebted to Mr. Diwan 
for the opportunity of emphasizing its importance. 


11. Ibid. at 10. 
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FIGURE 14. SLOPES AND DEFLECTIONS AS AFFECTED BY REDUCTION 
IN STEEL STRAIN DUE TO CONCRETE COVER 


| 
700 
L 
0.0006 & 0.0003 
|] mi 
2 
0.0008 0.0004 
0.2L 0.4L 0.4L 0.2L 
813-20 


FIGURE 


1S. 


FIXED END BEAM 
UNIFORM LOAD 
EQUAL POS. & NEG. STEEL Ph 
SPAN LENGTH=L 
f,» = 40,000 


0.04L 0.08 L 0.12L 0.164 0.20L 


min 


STEEL RATIO VS. MINIMUM LENGTH OF DIS~- 
TRIBUTION (Amin) FOR EQUAL POSITIVE AND 
NEGATIVE STEEL RATIOS. 


0.07 
= 
0.05 
0.04 
0.03 
Pp 
0.01 
0 
813-21 


at 


Discussion of 
“LATERAL BENDING OF SUSPENSION BRIDGES” 


by Cevdet A. Erzen 
(Proc. Paper 663) 


A. SELBERG.*—Since 1932,when Moisseiff and Lienhard(4) published their 
paper on the lateral bending of suspension bridges, this problem has been fre- 
quently handled in the literature. | -10 However, already in 1906 an approxi- 
mate solution was given by J. Melan(1), and description of a more exact in- 
vestigation was given on the same place. 

Investigations of the same problem were also given in 1919(2) and 1921.(3) 

The very interesting paper by C. A. Erzen gives a method which may be 
used on all possible types of suspension bridges. The method is, however, 
laborious, and takes a lot of unnecessary time for the designer. 

To get a quick and accurate solution of this problem, the best method is to 
express the horizontal component of tension in suspender ropes T as a 
series.(6) For instance 
x x 
1 1 


For symetrical spans Ag; Ag etc. will be zero. 

For smaller suspension bridges the terms A, or Aj ; Ag will be sufficient 
for all practical purposes, and for greater spans Aj; Ag; Ag usually will be 
sufficient. Calculating the example used by Moiseiff and Lienhard(4) with the 
terms Aj; A3 gives a difference in moments less than 2%, with 3 terms less 
than 1%. If T is expressed by a series u, M and v are found from equations 
(1), (3) and (5) by integration. 

The simplest method for finding the coefficients in the series for T is to 
use the equation (4). The equation shall be satisfied for so many values of x 
as there are coefficients in the series. 


+ Ag ~+- - - 


T = Ay sin7 i 


+ Ag sin 27 
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CEVDET A. ERZEN.*—The writer is much grateful to Dr. Selberg for his 
discussion of the problem of bending of suspension bridges under the effect of 
wind forces. The extensive literature on the problem as pointed out by Dr. 
Selberg is noteworthy. 

It is interesting to see the method by which this problem is solved by using 
a finite trigonometric series. The method, known as the Collocation Method, 
has been used quite extensively in the solution of differential equations. How- 
ever, it is necessary that each term of the series chosen for the purpose must 
approximate in shape the variation of the function that it represents. Other- 
wise convergence is not easily attained. 

The trigonometric series recommended in the discussion is apparently one 
that gives very rapid convergence if, with only two terms, one can obtain re- 
sults within two percent. Integration of Eqs. (1) and (3) with T expressed by a 
trigonometric series is a simple matter if py is constant throughout the span. 
However, it may take a considerable time to integrate p; if it is a function of 
x as given in the example. In this case, the writer believes that using a finite 
trigonometric series approach is not a simpler way to solve the problem, 


* Asst. Prof. of Civ. Eng., Cornell Univ., Ithaca, N. Y. 
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